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ABSTRACT 

We explore the possible connection between the open cluster IC 2391 and the unbound Ar- 
gus association identified by the SACY survey. In addition to common kinematics and ages 
between these two systems, here we explore their chemical abundance patterns to confirm if 
the two substructures shared a common origin. We carry out a homogenous high-resolution 
elemental abundance study of eight confirmed members of IC 239 1 as well as six members 
of the Argus association using U VES spectra. We derive spectroscopic stellar parameters and 
abundances for Fe, Na, Mg, Al, Si, Ca, Ti, Cr, Ni and Ba. 

All stars in the open cluster and Argus association were found to share similar abundances 
with the scatter well within the uncertainties, where [Fe/H] = -0.04 ± 0.03 for cluster stars 
and [Fe/H] = -0.06 + 0.05 for Argus stars. Eff'ects of over-ionisation/excitation were seen 
for stars cooler than roughly 5200 K as previously noted in the literature. Also, enhanced Ba 
abundances of around 0.6 dex were observed in both systems. The common ages, kinematics 
and chemical abundances strongly support that the Argus association stars originated from 
the open cluster IC 239 1 . Simple modeUng of this system find this dissolution to be consistent 
with two-body interactions. 

Key words: (Galaxy): open clusters and associations: general - (Galaxy): open clusters and 
associations: individual: IC 2391 and Argus association - stars: abundances 



1 INTRODUCTION 

The vast majority of stars that we observe in any given galaxy are 
field stars, a complex mixture of stellar populations of different 
ages and metallic! ties. The present view is that star formation 
taking place in molecular clouds gives rise to groups of stars 
- stellar associations and open clusters. These groups of stars 
must dissolve with time to connect star formation with the field 
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population. 

Based on ne ar-infrared studies of embedded clusters, 
iLada & Lad3 ( l2003h suggest that less than about 10% of the clus- 
ters formed in molecular clouds survive longer than 10 Myr. They 
propose that this "infant mortality" of clusters is due to feedback 
of massive stars that expels the gas from the star forming region 
leaving the stars behind with super-virial velocities. For compact 
systems of a few hundred stars the dynamical time-scale is only a 
few Myrs. This means that cluster expansion and star loss by two- 
body rel axation is also a relevant process at th ese ages for such 
systems jMoeckel etaTlbOljjGieles et al.ll2012h . On longer time- 
scales (few 100 Myrs) various processes such as the Galactic tidal 
field and interactions with Giant Molecular Clouds (GMCs) be- 
come important. The paucity of open clusters in the solar neigh- 
bourhood with ages between 0.5 and a few Gyrs was already noted 
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in the late 50s toortlll958h . lSpitzeil figSSh analytically derived the 
cluster disruption times due to encounters with GMCs, showing 
that such phenomena would account for the scarcity of Galactic 
clusters older than 10' years. A more detailed treatment of the prob- 
lem taking into acco unt other disruptive agen ts has been carried out 
in th e literature (e.g.lLamers & Gielesll2006l) . 

[Torres et all ( |2006|) reported the results of a high-resolution 
optical spectroscopic survey to search for associations containing 
young stars (SACY) among optical counterparts of ROSAT All- 
Sky X-ray s ources in the southern hemisphere. Us ing the method 
described in lTorres et al. I (l2006l) .l Torres et al. identified nine 

new young associations, namely, Chamaleontis (ChA), TWHydrae 
(TWA), 13 Pictoris, Octants (OctA), Tucana-Horologium (THA), 
Columba (ColA), Carina (CarA), Argus (ArgA), and AB Doradus 
(ABDA). The age span derived based on the isochr one fitting and 
the L i abundance analysis ranges from 6-70 Myr ( Ida Silva et alj 
l2009l) . 

Of relevance to the context of cluster disruption discussed 
above, are two young associations identified in the SACY survey 
which appear to be associated with known open clusters. These 
are the Argus association and the open cl uster IC 2391 and the 
e Cha Association and the rj Cha cluster jMamaiek et"ani200(]| : 
[Torres et al.l |2008|) . In this paper we investigate in detail the 
possible connection between the Argus Association and the young 
open cluster IC 2391. 

In particular, we can test for a common origin of the associ- 
ation and cluster by examining the stellar elemental abundances. 
If the abundance analysis reveals high chemical homogeneity 
among the association stars comparable to the open cluster 
abundance patterns, it supports a common origin scenario where 
the Argus association stars were once bound members of open 
cluster IC 2391. This assumes that the proto-cluster gas cloud 
was sufficiently mixed hence all stars bom from the same site 
should share the same elemental abundance pattern. This concept 
of elemental abundance matching is often referred to as chemical 
tagging (Freeman & Bland-Hawthom 2002s). Conversely if the 
abundances of the Argus association stars are found to be either 
heterogeneous, matchin g the general disk s tellar population (e.g. 
the Hercules stream, see lBensbv et al. I( l2007h ) or they do not match 
the abundance patterns of IC 2391, then it is more likely that these 
are two separate systems even if they share the same kinematics. 

Given that disk dynamics make the kinematical information 
somewhat unreliable to explore the history of these stellar struc- 
tures, for the first time we carry out a high resolution elemental 
abundance analysis of the Argus association in a homogeneous 
maimer with stars in IC 2391 in order to carry out chemical tag- 
ging. This will also be the first attempt at chemical tagging a young 
stellar system. All published examples of chemical tagging in the 
disk are for stellar systems older than the Hyades dBubar & Kind 
l20I(]| : lDe Silva et alfcoOTbl) . 

In Section[2l we briefly describe the SACY survey and present 
space velocities and spatial locations of the stars in the two sys- 
tems. In Section[3]we present a detailed elemental abundance study 
of stars in both the Argus association and open cluster IC 2391. 
In Section |4] we discuss our results in terms of chemical tagging 
the two stellar systems. We conclude with a discussion of possible 
disruption mechanisms given the evidence gathered in the previous 
sections. 



2 IDENTIFYING THE ARGUS ASSOCIATION IN THE 
SACY SURVEY 

2.1 The convergence method 

Usually we think about an association as being a group of stars 
appearing concentrated in a small volume in space sharing some 
common properties such as age, chemical composition, distance 
and kinematics. However, if such a group is close enough to the 
Sun, its members will appear to cover a large extent in the sky (as 
an example, Orion at 50 pc would cover almost the whole sky). 
Thus, to find a group, projected spatial concentrations (i.e., in 
terms of right ascension and declination only) and proper motions 
may not be enough. A better criterion is to look for objects sharing 
similar heliocentric space motions (UVW) all around the sky 
(where we use U positive towards the Galactic center, V positive in 
the direction of Galactic rotation, and W positive in the direction 
of the North Galactic Pole). 

[Torres et al. I ( l2006l) describe in some detail the convergence 
method developed to search for members of an association. 
This method examines the stars in the hexa-dimensional space, 
UVWXYZ, as defined by the space motions relative to the Sun and 
the physical space coordinates centered on the Sun (XYZ, in the 
same directions as UVW). We represent with m,,, M,, and M,,,,^ 
the apparent visual magnitude, the resultant absolute magnitude 
with the distance obtained from the convergence method, and the 
absolute magnitude given by the adopted isochrone for the (V - /)c 
stellar color; and yu„, yu,; and V,- are the proper motions and the ra- 
dial velocity. In brief, if there is no reliable trigonometric distance 
availablfl the convergence method finds the distance (d) for each 
star in the sample that minimizes the F value of Equation [T] The 
first term is a photometric distance modulus and the second is a 
kinematical one. The method needs, as input, an assumed age and 
initial velocity values {Uq, Vq, Wq) for the proposed association, 
and a cutoff value for F above which stars should be considered 
spurious. This cutoff value varies for each association but usually 
we begin with 3.5 (this approximately means 0.7 magnitudes for 
the distance modulus and 3 km s"' for the velocity modulus). 

The method is iterative, and for each iteration a list of stars 
with new (t/o, Vo, Wo) is obtained. The process ends when the list 
of stars and the velocities (t/o, Vo, Wo) do not change significantly. 

F(mi,,yU„,yU,s, Vr\d) = 

[p X (M., - M„,,,„)2 + (U- Uof + (V- Vo? + {W- Wcf-VHl) 

where p is a constant weighting the importance of the evolutionary 
distance with respect to the kinematic distance. Usually, if there are 
several stars with trigonometric parallaxes among the candidates, 
we use p=0. But if there are no stars with reliable parallaxes, as is 
the case for the IC 2391, the method may not converge, so we used 
p=20 for this case. 

We also define a membership probability to verify the com- 
pactness of the association in 7 dimensions, namely UVWXYZ and 
(M,,-M,, , ,„). Normalized values of these variables are computed by 
subtracting their mean values and dividing by their dispersions. A 



We consider the trigonometric parallaxes as unreliable if they have errors 
larger than 2 mas. 
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Figure 1. Combinations of the UVWXYZ-space derived based on the convergence method for the Argus Association stars (red squares) , the open cluster 
IC 2391 (blue triangles) and when the two systems are considered as a single association (gray circles). Well defined clustering can be seen in both kinematical 
and spatial coordinates, with no significant difference when considering the two stellar system independently or together. 



quantity K is obtained as the quadratic sum of these dimensionless 
variables and the probability is computed as 

P=l- {ermy (2) 

where erf is the error function. 

Finally, we also take into account the Li content of the candi- 
date stars to see i f it is compatible with the Li depletion for its age 
( lNeuhausej|l997h . 
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Table 1. Revised list of high probability members for the Argus association 



Name 


a 


6 


Ma 






Vsini 




Vj 


V-Ic 


SpT 


El, 




(V, V. W) 


(X, Y, Z) 




Prob 


notes 




2000 


2000 


mas yr ' 


km 




days 








mA 


mas 


km s ' 






pc 








(*)CD-29 2360 


05 34 59.2 


-29 54 04 


1 7 1 


3 J. I 


26.0 




3.519 


10.53 


1.11 


K3Ve 


180 


14.8 


-22.7-16.6 


-5.2 


-34.9 


-47.9 


-32.4 


0.72 


K 


(*)AP Col 


06 04 52.2 


-34 33 36 


27 3 


340 9 


22.4 


n 




12.96 


3.36 


M4Ve 


280 


119.2 


-22.0-13.6 


-4.4 


-3.7 


-6.7 


-3.4 


0.90 


1 


CD-56 1438 


06 11 53.0 


-56 19 05 


-2 5 


JO. J 


14 


130 


0.418848 


11.32 


0.93 


KOV 


230 


8.7 


-21.9-11.2 


-5.1 


-9.1 


-101.3 


-53.3 


0.68 


K 


CD-28 3434 


06 49 45.4 


-28 59 17 


1 1 1 
-I I. J 


zU.o 


26.8 


6.4 


3.823 


10.62 


0.85 


G7V 


230 


9.8 


-23.1-16.3 


-7.0 


-51.2 


-85.3 


-23.1 


0.45 


K 


CD-42 2906 


07 01 53.4 


-42 27 56 


1 ^ 8 




23.6 


10.8 


3.96 


10.64 


1.02 


KIV 


275 


10.7 


-22.8-16.8 


-6.3 


-26.5 


-85.8 


-26.0 


0.75 


K 


CD-48 2972 


07 28 22.0 


-49 08 38 


-zo. / 


AA A 


21.1 


52.0 


1.0373 


9.85 


0.86 


G8V 


250 


12.6 


-22.3-16.6 


-7.3 


-11.9 


-76.0 


-20.0 


0.88 


K 


(*)HD 61005 


07 35 47.5 


-32 12 14 


-J J. / 


7A A 


22.5 


8.2 


5.04 


8.22 


0.80 


G8V 


171 


28.3 


-23.0-14.1 


-4.4 


-14.1 


-32.2 


-3.5 


0.73 


H,2 


CD-48 3199 


07 47 26.0 


-49 02 51 


-22 4 


36 1 


18.5 


24.7 


2.192 


10.56 


0.86 


G7V 


230 


10.0 


-22.2-15.1 


-5.0 


-12.9 


-96.7 


-20.2 


1.00 


K 


CD-43 3604 


07 48 49.8 


-43 27 06 


-ZO.o 


Jy.v 


21.4 


40.4 


0.4706 


11.13 


1.40 


K4Ve 


320 


12.2 


-22.6-16.4 


-4.8 


-17.5 


-79.1 


-12.6 


0.82 


K 


TYC 8561-0970-1 


07 53 55.5 


-57 10 07 


- 1 O.J 


7^^ Q 
ZD.y 


15.8 


5.0 




11.50 




KOV 


210 


7.0 


-21.9-14.7 


-6.5 


0.6 


-139.0 


-36.4 


0.98 


T 


HD 67945 


08 09 38.6 


-20 13 50 


^7 R 
-J 1 .o 


7n A 

ZU.4 


22 


120 




8.08 




FOV 





14.3 


-22.1-13.4 


-4.2 


-34.7 


-60.0 


8.5 


0.84 


TSB2? 


CD-58 2194 


08 39 11.6 


-58 34 28 


^7 7 


A 

J J. 4 


15.6 


85.0 


1.240 


10.18 


0.83 


G5V 


270 


9.8 


-21.9-16.7 


-5.4 


8.3 


-100.3 


-18.3 


0.98 


K 


CD-57 2315 


08 50 08.1 


-57 45 59 




^ 1 n 


11.7 


52.0 




10.21 




K2Ve 


308 


9.3 


-22.1 -13.0 


-5.5 


9.2 


-105.8 


-16.2 


0.90 


T 


TYC 8594-0058-1 


09 02 03.9 


-58 08 50 


7A 7 


If^ 7 
ZD.Z 


12.1 


34.1 




11.30 


0.83 


G8V 


300 


7.2 


-22.2-14.4 


-2.4 


15.2 


-137.7 


-18.7 


0.95 




CD-62 1197 


09 13 30.3 


-62 59 09 


'XA 1 


77 A 
Z / .4 


12.7 


84.0 


0.54841 


10.46 


0.93 


KOV(e) 


280 


8.4 


-21.7-16.3 


-5.8 


22.1 


-114.9 


-20.4 


0.91 


K 


TYC 7695-0335-1 


09 28 54.1 


-41 01 19 


7R 7 


^J.\J 


14 


120 


0.3917 


11.65 


1.10 


K3V 


300 


6.7 


-22.0-13.6 


-5.3 


-7.9 


-147.3 


18.8 


0.87 


K 


(*)BD-20 2977 


09 39 51.4 


-21 34 17 


- JU.O 


1 8 


18.3 


10.1 




10.22 




G9V 


260 


11.3 


-22.3-16.5 


-4.7 


-21.9 


-79.0 


34.4 


0.58 


T 


TYC 9217-0641-1 


09 42 47.4 


-72 39 50 


'^n 1 

- JU. i 




6.7 


23.4 


2.3026 


12.30 


1.11 


KIV 


240 


6.5 


-22.1 -13.6 


-6.2 


50.4 


-139.6 


-39.1 


0.88 


K 


CD-39 5833 


09 47 19.9 


-40 03 10 




ID. J 


15.0 


10.5 




10.89 




KOV 


260 


9.0 


-22.2-15.7 


-4.6 


-2.1 


-109.5 


20.0 


0.92 


T 


HD85151A 


09 48 43.2 


-44 54 08 


A 


'X'X ^ 
J J.J 


14.3 






9.61 




G7V 


220 


15.3 


-22.3-15.8 


-3.9 


2.5 


-64.8 


7.7 


0.92 


T 


HD85151B 


09 48 43.4 


-44 54 09 


A 


J J. J 


13.4 






10.21 




G9V 


250 


15.3 


-22.3-14.9 


-4.0 


2.5 


-64.8 


7.7 


0.84 


T 


CD-65 817 


09 49 09.0 


-65 40 21 


?7 ? 
-JZ.J 


/O 


7.3 


19.2 


2.7388 


10.17 


0.75 


G5V 


200 


7.1 


-22.4-13.1 


-4.8 


37.5 


-134.5 


-22.5 


0.97 


TK,D2.0" 


HD 309851 


09 55 58.3 


-67 21 22 


AO A 
-4Z.4 


7< 7 
ZJ. / 


6.6 


19.8 


1.8164 


9.90 


0.68 


GIV 


170 


9.2 


-22.5-13.2 


-4.5 


31.8 


-102.6 


-19.0 


0.98 


K 


HD 310316 


10 49 56.1 


-69 51 22 


-43.5 


13.5 


4.7 


16.3 


3.625 


10.07 


0.78 


G8V 


224 


8.2 


-22.4-13.6 


-5.2 


46.5 


-110.7 


-19.9 


0.97 


K,D0.6" 


CP-69 1432 


1053 51.5 


-70 02 16 


-32.6 


12.2 


3.6 


55.0 


1.0305 


10.66 


0.71 


G2V 


195 


6.3 


-22.8-13.1 


-3.3 


61.8 


-144.6 


-26.1 


0.64 


K 


(*)CD-42 7422 


12 06 32.9 


-42 47 51 


-50.3 


-2.5 


3.3 


28.3 


1.984 


10.66 


0.92 


KOV 


260 


8.9 


-21.6-15.4 


-4.8 


43.4 


-96.5 


37.1 


0.50 


K 


CD-74 673 


12 20 34.4 


-75 39 29 


-109.4 


4.2 


2.7 


7.3 


3.477 


10.72 


1.17 


K3Ve 


230 


19.5 


-21.7-15.5 


-2.8 


25.8 


-42.9 


-11.5 


0.75 


K,3,4 


CD-75 652 


13 49 12.9 


-75 49 48 


-62.2 


-31.5 


-0.9 


20.1 


2.2743 


9.67 


0.76 


GIV 


200 


12.3 


-22.3-13.6 


-6.0 


47.0 


-63.4 


-18.8 


0.93 


K 


HD 129496 


14 46 21.4 


-67 46 16 


-47.6 


-39.7 


-4.8 


88.0 




8.78 


0.60 


F7V 


150 


11.1 


-22.5-13.6 


-5.9 


61.6 


-64.9 


-11.5 


0.79 




NY Aps 


15 12 23.4 


-75 15 16 


-73.9 


-73.1 


3.4 


10.8 


4.084 


9.42 


0.86 


G9V 


182 


19.9 


-21.2-12.6 


-4.3 


32.3 


-36.3 


-12.9 


0.62 


K,H 


(*)HD 145689 


16 17 05.4 


-67 56 29 


-50.0 


-84.0 


-9 


106.4 




5.95 


0.18 


A6V 




19.2 


-22.7-11.7 


-3.7 


39.2 


-32.4 


-11.2 


0.68 


H,5.6,7 


CD-52 9381 


20 07 23.8 


-51 47 27 


85.9 


-143.9 


-13.3 


42.0 


0.8368 


10.59 


1.52 


K6Ve 


60 


33.6 


-22.4-14.8 


-4.2 


24.4 


-5.7 


-16.0 


0.80 


K 



'Data in italics are fr om the literature as follo ws: H = data from Hip parcos; T = data from TY CHO; K = data fromlKiragal |201^: (niRiedel et al.l i201lh ; (2) iDesidera et al.l j201 lah ; (3) SBl, P<„-(,=613.9d 
iGuentheret ai] j2007l) ; (4) ICovino et"ai] l ll997h : (5) Brown d warf at 6.7"lHiielamo et alj j2010l) ; (6) IZuckerman et al.l \20\ ih : (7) lDi'az et al.l j201lh . 
Entries marked (*) are new inclusions previously not listed bv lToiTes et al.l j2008l) . 



Table 2. Revised list of high probability members proposed for the open cluster IC 2391 



Name 


a 


S 






Vr 


Vsini 


P/;/;or 


Vj 


B-Vor 


SpT 


El, 




(V 


, V, W) 




(X,F,Z) 


Prob 


notes 




2000 


2000 


mas 


yr"' 


km 


S-' 


days 




V-Ic 




mA 


mas 


km s ' 






pc 








(*)PMM 5314 


08 28 34.6 


-52 37 04 


-24.1 


22.7 


13.4 


63.0 




10.40 


0.45 






7.1 


71 7 


1 9 9 
-i Z.Z 


-6.3 


-2.5 


-138.5 


-19.6 


0.92 




PMM 7422 


08 28 45.6 


-52 05 27 


-22.6 


24.6 


14.8 


33.0 


1.515 


10.49 


0.69 


G6 


233 


7.2 


-22 1 


1 '^ 7 


-4.8 


-3.5 


-137.1 


-18.6 


1.00 


K 


PMM 7956 


08 29 51.9 


-51 40 40 


-23.0 


18.8 


15.3 


13.8 




10.62 


0.98 


e 


289 


6.5 




1 ^ 8 
-i J.o 


-7.7 


-4.6 


-152.7 


-19.7 


0.71 




PMM 6974 


08 34 18.1 


-52 15 58 


-22.6 


23.0 


15.2 


5.2 


7.80 


12.26 


1.04 




74 


6.9 


90 1 




-5.1 


-2.1 


-143.5 


-18.0 


1.00 


1 


PMM 4280 


08 34 20.5 


-52 50 05 


-21.0 


22.5 


12.3 


16.0 




10.34 


0.67 


G5 


151 


6.6 


-22 1 


-117 


-4.3 


-0.9 


-151.0 


-19.8 


0.90 




PMM 6978 


08 35 01.2 


-52 14 01 


-23.3 


22.3 


15.2 


7.1 


5.1 


12.07 


1.02 




179 


6.9 


-zz.u 


1 A 


-5.7 


-2.0 


-143.6 


-17.7 


1.00 


1 


PMM 2456 


08 35 43.7 


-53 21 20 


-23.6 


24.3 


14.8 


48.0 




12.20 


0.92 


K3e 


301 


7.2 


99 1 
-ZZ. i 


1 A 


-5.0 


0.5 


-138.0 


-18.4 


1.00 




PMM 351 


08 36 24.2 


-54 01 06 


-24.4 


20.7 


15.4 


90.0 


1.923 


10.18 


0.57 


GO 


90 


6.8 


91 fi 
-Zi .o 


1 A S 
-i'+.o 


-7.3 


2.0 


-146.2 


-20.3 


0.92 


1 


PMM 3359 


08 36 55.0 


-53 08 34 


-23.1 


23.6 


14.7 


8.2 
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2.2 Kinematics and ages of IC 2391 and the Argus 
Association 

The Argus association was easily discovered in tiie SACY survey 
due to its speci al U velocity of U = -22.0km s"' as illustrated 
in Figure 2 of [Torres et alj ( 1200 8l) . The evidences presented in 
[Torres et"an l l2008h stronglv suggested that the Argus association 
and IC 2391 have a common origin. We have updated the SACY 
catalog in order to be able to test if memb ers proposed by 
iDesidera et"ZI ( l2011ah . IZuckerman et al.l ( l201lh . and iRiedel et all 
1 201 ih could be converged in to the IC2391 and Argus association 
member list. Newly obtained data presented in this paper were also 
used. 



The c onvergence was do ne using the new reduction of 
Hipparcos (van Leeuwen '2007*), the new proper motion catalog 
UCAC4 (iZacharias et al. 2012) and the AS AS photometric survey 
(lKiragjl2012h . for V and V-I colors, for stars not observed by us. 
We did not change the published proper motions for stars in the 
Hipparcos catalogue, and for double stars in UCAC4 with sepa- 
rations between about 1" to 10" because the image de-blending 
proce ss in UCAC4 seems unreliable. As we did in dTorres et al.l 
l2008h . we correct the photometric values for duplicity where 
possible. Note for the pair HD 85151A/B, there are no individual 
proper motions in Tycho-2. Therefore we used the primary proper 
motions for both stars. Our spectroscopic observations and the 
visual binary measurements from the Washington Double Star 
catalogue suggest they must be a physical pair. 



The data for IC 2391 members are mainly from lPlatais et al.l 

( l2007h but usi ng UCAC4 proper moti ons, except for doub le stars 
where we use IPlatais et"al] ( l2007h and iTorres et al ] (l2008h proper 
motions due to image de-blending issues with UCAC4 as explained 
above. There are 41 cluster stars with all kinematical data. Except 
for PMM 4413 whose systemic velocity is known, all other 20 
known spectroscopic binaries were excluded from the study. 

In Table [T] and Table [2] we list the high probability members 
proposed for the Argus association stars and the open cluster 
IC 2391 respectively, obtained by applying the convergence 
method for both association and cluster stars together. We present a 
total of 65 high probability members which comprise of 32 Argus 
association stars and 33 IC 2391 members. We accept probabilities 
down to about 50 percent as the stellar system is formed by the 
concentrated IC 2391 cluster and the widely scattered association 
stars around it, thus forming a spatial distribution that is far 
from gaussian. When we apply the convergence criteria on the 
association stars alone (without including IC 2391 members), then 
the probabilities of the association members increase significantly 
with the lowest membership probability being about 70 percent.) 

In order to test the true kinematical connection between 
IC 2391 and the Argus association stars, we also applied our 
convergence test to the two groups independently. The results 
are presented in Table |3] and in Figure [T] 1 The convergence 
method resulted in a near perfect match between the kinematics 
of IC 2391 stars a nd those of the Argus association stars from 
[Torres et"al1 (l2003h . The velocities of the two samples agree well 
within the errors. To give a better idea of the variations in the final 
UVWXYZ values when the Argus association stars are considered 
together with the IC 2391 members, we have over-plotted the 
results obtained for the Argus field stars only (red squares), and 




Figure 2. The HR diagram of the members proposed for the Argus Asso- 
ciation. Blue triangles are the ArgA field members whereas the IC 2391 
members are rep resented by red triangles. The overplotted isochrones 
iSiess et alj200ol) are the ones for 20 (dot-dashed) , 30 (solid) and 40 Myr 
(dashed). 



the IC 2391 members only (blue triangles), with those obtained 
when the two groups are considered as one single association (gray 
circles). No significant difference is found. 

Also as shown in Table |3] no significant difference is found 
in the ages of the two groups. Figure |2] shows that the position of 
the Argus association in the HR diagram is consistent with a unique 
isochrone of 30 Myr +10 Myr, which is comparable to our obtained 
age of about 26 Myr. Our ages are from third degree polynomial 
isochome curves dTorres et al.lbOOSi) based on some well observed 
star groups and compared partially with lSiess et al. I (l2000h . As they 
are not strictly calibrated, we must take these heuristic ages with 
caution, but of course the relative ages are more reliable. These 
results imply that bo th samples have the same age. Additionally, 
Ida Silva et atl ( l2009h computed Li abundances for all SACY asso- 
ciations and showed that the Lithium Depletion Pattern for the Ar- 
gus association (defined by fitting the Li abundance as a function 
of effective temperature), fits well in a time sequence of the SACY 
associations with no difference between the Argus association stars 
and the IC 2391 members. 



2.3 Conunents on excluded stars 

Table [T] and Tabl e |2] are not the same as those presented in 
iTorres et al. I (l2008l) . as we used updated data and applied a new 
statistical model. Below we comment on the excluded stars. 

For th e Argus asso c iation sample, 3 stars were excluded from 
the list of ITorres et ail ( |2008|) . Star BW Phe was eliminated on 
kinematical grounds, star CD-49 1902 was eliminated for having 
too large value in Z and star HD 84075 has a U velocity well 
outside of the Argus association range. 

For the IC 239 1 sample, 3 stars were excluded from the list of 
ITorres et al.l ( |2008|) . PMM 1560 and PMM 1820 were eliminated 
due to low membership probabilities based on kinematics. For 
PMM 3695 the proper motions are significantly distinct between 
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Table 3. Convergence method results for the Argus field stars and the members of open cluster IC 2391 . 



Sample Distance [mas] U [km/s] V [km/s] W [km/s] N. members Age [Myr] 

Argus only (field) 15.5 ± 19.9 -22.6 ± 0.4 -14.6 ±1.5 -5.0 ± 1.2 32 26 

IC2391 6.9 ±0.4 -21.7 ±0.6 -14.2 ±1.0 -5.1 ± 1.2 33 26 

Argus &IC 2391 11.2 ±14.5 -22.1 ±0.4 -14.4 ±1.3 -5.0 ± 1.2 65 26 



IPlatais et alj ( l2007h and UCAC4, for unclear reasons. 



3 ELEMENTAL ABUNDANCES 

3.1 High resolution observations 

In order to further explore the connection between the open cluster 
IC 2391 and the Argus association, we examine their elemental 
abundances. High resolution and high signal-to-noise spectra of 
members of the Argus association and IC 2391 were observed 
using VLT-UVES in the framework of program ID 082.C-0218 (PI 
Melo). Observations were undertaken in service mode using the 
UVES DIC#1 (390-l-580nm) standard setting with a O.Barcsec slit 
to achieve a spectral resolving power of R=60,000. The typical 
S/N ratio was ~100 per pixel at 600 nm. The data were reduced 
with the latest UVES ESO-MIDAS pipeline. The resulting spectra 
were normalised using the continuum task in the IRAlQ package. 

We carry out a detailed elemental abundance analysis of six 
Argus association stars and eight IC 2391 stars with vsini < 10 
kms"'. The other stars showed broadened spectral lines which 
meant there was greater blending between the lines, making indi- 
vidual line analysis more uncertain. Table|4]lists the sample of stars 
analysed in this study together with their measured properties. 

3.2 Abundance analysis 

The elemental abundances were derived based on equivalent width 
(EW) measurements and spectral synthesis, ma king use of the 
latest version of the MOOG code l lSneden|[T97|)^ The EWs were 
measured using the automated ARES code Isousa et alj l2007h 
with frequent manual checking of the EWs. Ba abundances were 
derived from the 5853 A line (which is known to be strong, 
isolated, un-blended an d not affected by NUTE effects, e.g., 
iMashonkina et alj ( |2007|) ) using the dri ver synth in MOOG. \V e 
retrieved hyperfine structure data from iMcWilham eTaTI (Il995h . 
adopting an isotopic solar mixture of 81% ('^*Ba -l-'^*Ba -l-'^^Ba) 
and 19% ("^Ba -l-"^Ba). An example of the synthesis of the Ba ii 
line at 5853 A is provided in Figure |3]for one of our sample stars 
(PMM 3359). 

Interpol ated Kurucz mode l atmospheres based on the AT- 
LAS9 code JCastelh etal]|l997l) with no convective overshooting 
were used throughout this study. We carried out the abundance 
analysis relative to the solar spectrum, taken with the same 
instrument and the same resolution of our sample stars. The full 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 



lines list, atomic data adopted and measured equivalent widths of a 
sample star (PMM 3359) can be found in Appendix IaI 



We derive the stellar parameters based on spectroscopy. Abun- 
dances for all Fe i and ii lines were computed from the measured 
EWs. Elfective temperature (Jeff) was derived by requiring exci- 
tation equilibrium of the Fe i lines. Microturbulence was derived 
from the condition that abundances from Fe i lines show no 
trend with EW. Suiface gravity (log g) was derived via ionisation 
equilibrium, i.e. requiring the abundances from Fe i lines to equal 
those from Fe ii lines. For stars cooler than 5500 K, we adopted 
a logg value of 4.5, which is consistent with the cluster age. We 
did not attempt to optimize logg via ionisation equilibrium for 
the cool stars becau se of the known over-ionisation effects (e.g., 
ISchuler et al.l ( l201Clh . see Section |34l for further discussion). The 
final adopted stellar parameters and the elemental abundances are 
given in Table |4] 



3.3 Error budget 

Chemical abundances are mainly affected by two kinds of un- 
certainties, namely (i) errors due to the EW measurements or to 
the best-fit determination (in the case of spectral syntheses), and 
(ii) errors due to the stellar parameters (Teff, logg, and ^). The 
impact of the uncertainties on the atomic parameters, i.e. log gf, 
should be instead almost negligible since our analysis is strictly 
differential with respect to the Sun (see Section 13.2) . Concerning 
the EW analysis, random errors in the [Fe/H] ratios due to EW 
measurements are well represented by the standard deviation (rms) 
from the mean abundance based on the whole set of lines, while 
the ones for the [X/Fe] ratios were computed by quadratically 
adding the rms for [Fe/H] and rms for [X/H]. Errors due to stellar 
parameters were estimated by varying one parameter at a time, and 
checking the corresponding variation in the resulting abundance. 
We adopted variations of ±50 K in Ted and +0.15 kms"' in 
because larger changes in those quantities would have introduced 
a significant trend of log n(Fe) vs the excitation potentials and the 
line strength, respectively. For stars for which we could optimise 
the surface gravities, the uncertainties in log g were estimated by 
varying this quantity until the difference between log ;!(Fei) and 
log ;!(Feii) is larger than 0.1 dex, i.e. the ionisation equilibrium 
condition is no longer satisfied. The typical error in log g was from 
0.1 to 0.15. The final enws are then obtained by quadratically 
summing errors due to the EW measurements and the ones due to 
stellar parameters (see Table |4). 

Focusing on abundances derived from spectral syntheses, we 
have errors due to the best-fit determination of typically 0. 1 dex, 
reflecting uncertainties due to the continuum placement and to the 
strong nature of the Ba features; similarly, since we are dealing with 
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Figure 3. Example of the Ba synthesis for star PMM 3359, with solid red line shows the best fit. 



a very strong transition line, the errors due to stellar parameters are 
dominated by the microturbulence values (A[Ba/Fe]=0.05-0.07 dex 
for ^ changes of 0.15 kms"'), while the effective temperatures and 
gravities have a minor impact, i.e. within 0.03 dex. 



3.4 Abundance results 

Final abundances are reported in Table |4] where stellar parameters 
(Teff, log^, ^) and [X/Fe] ratios along with the corresponding 
errors are provided for all our sample stars. In Figure 2]we show 
the [X/Fe] ratios as a function of the effective temperatures for 
both IC 2391 (filled circles) and the Argus association (triangles). 
The first evidence c oming out from these plot s is that we confirm 
previous findings bylSchiiler et al.l(l2003ll2004l):IShen et al.l(l2005h; 



Ip'Orazi & Randiclj ( EoogI) : ISchuler et al.l ilOld) : ISiazzo et alj 
1201 l|): over-ionisation/excitation effects are at work for stars 
cooler than Teff<5200 K. The photospheric abundances of Ca, 
Na, and Ti (from the Ti i lines) show the most dramatic effect of 



this phenomenon, being about 0.2-0.3 dex lower in the cooler 
stars; similarly, Al and Cr (from the Cr i lines) seem to show 
the same effect, although to a lesser extent. Because we could 
also measure abundances for Cr and Ti from their first-ionisation 
stage transitions, we can confirm that we are observing an over- 
ionisation phenomenon: as can be seen in Table [4] abundances 
from Ti ii and Cr ii lines result in larger values for the cooler 
stars, where the [X/Fe] ratios are enhanced up to 0.4 dex (see the 
[Cr ii/Fe] in the most extreme cases, e.g., PMM 4902 in IC 2391 
or TYC 8561-0907-1 in Argus association). This fact reflects 
a pumping of the electrons to the ionised stage because of the 
strong U V flux coming from the hot c hromospheres (we refer the 
reader to ISchuler et all (l2003l |2010|) : iD'Orazi & RandichI fcOOSl) 
for a more detailed discussion on this topic). Si deserves special 
attention since we found an over-excitation effect (all the employed 
lines have high excitation potentials, i.e. 4.92 < ^ <5.98 eV) 
acting on the cooler stars, that was also pre viously detected by 
ISchuler et aO ( |2003|) in M34, but not by IP'Orazi & RandichI 
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( l2009h . Over-ionisation/e xcitation effects are a l so obs erved in the 
cool field stars (see e.g., lAUende Prieto et all ( 1200 Sh V Note that 
no temperature effects were observed for the Ni abundances in 
either the IC 2391 or the Argus association stars. Therefore the Ni 
abundance is not plotted in Figure |4] 

Due to the over-ionisation effects, the mean abundances and 
the associated rms error were computed considering all the stars 
for Fe I, Mg, Ni, and Ba, while we considered only the warmer 
stars (T^ff> 5200K) when averaging [X/Fe] ratios for Na, Al, Si, 
Ca, Ti and Cr (see Table |4j. We found average metallicity of 
[Fe/H]=-0.04±0.03 and [Fe/H]=-0.05+0.04 for IC 2391 and 
Argus association, respectively: the cluster and the association 
exhibit a solar metallicity, with no evidence of internal scatter (the 
rms is significantly smaller than the observational uncertainties). 
Moreover, all the other a (Mg, Si, Ca, and Ti), odd-Z (Na, Al), and 
iron-peak elements show solar ratios, pointing out that the cluster 
and the association are indistinguishable concerning their chemical 
composition (see Section |4] for a discussion on the scientific 
implication of our results). 

Ba is notably over-abundant, with a mean of 
[Ba/Fe]=0.62±0.07 (IC 2391) and [Ba/Fe]=0.53+0.08 (Ar- 
gus association). The slightly higher variation across the two 
stellar aggregates mainly reflect the larger uncertainties affecting 
Ba abundances. The Ba line feature is close to the saturation regime 
in the curve of growth, and as a consequence it is highly sensitive 
to the adopted microturbulence values. The extremely high Ba 
content that we obtained for our sample stars is not surprising and 
confirms previous measurements (see Section |4). We stress, in 
passing, that over-ionisation effects can be ruled out as a possible 
explanation given that, at the effective temperatures of our sample 
stars, Ba is almost totally ionised. Further the possibility of NLTE 
effects is minimal as the 5853 l ine used in this study is n ot known 
to be afi'ected by NLTE efl:ects jMashonkina et alj ( l2007h ). 
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Table 4. Stellai' parameters and abundances (see section |3!4l f or details on the computation of the average values) 
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1.3 


-0.03+0.12 




-0.25±0.15 


-0.06±0.09 


-0.08±0.06 


0.09 ±0.08 


-0.25±0.11 


-0.32±0.17 


O.I2±0.14 


-0.07 ±0.15 


0.41 ±0.23 


0.03 ±0.09 
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PMM 6974 
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4600 


4.5 
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-0.21+0.14 
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-0.09±0.04 
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-0.14±0.13 
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-0.31 ±0.1 3 


0.11±0.13 


-0.10±0.09 


0.22±0.18 
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13.9 
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1.2 


-0.05+0.10 
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-0. 1 1 ±0.09 
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-0.14±0.11 


-0.37±0.13 


0.21 ±0.11 


-0.19±0.11 


0.44±0.15 


0.01 ±0.09 
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3.5 Comparison with previous abundance studies 

The chemical composition of IC 2391 has been previously deter- 
min ed by two different studie s, namely Platais et al. 2007 (P07) 
and lD'Orazi&RandichI ( |2009|) . 

In P07, the authors presented iron abundances for a sample 
of four G-type members, obtaining an average metallicity of 
[Fe/H]=0.05±0.05, which is in very good agreement with our 
estimate. Three out of four P07 stars are in common with our 
study, i.e., PMM 3359, PMM 665, and PMM 4362; the differences 
in effective temperatures are ATeff(our-p07)=100 K, 99 K, and 
34 K, respectively, while gravities agree within 0.05 dex. Micro- 
turbulence values agree very well for the star PMM 3359 (A^=0.03 
km s"') and PMM 4362 (A^=0.1 km s"'), while P07 derived a 
slightly lower value for the star PMM 665, i.e., 1.25 km s"' to be 
compared with our estimate of 1.60 km s"'. 

iD'Orazi & RandiclJ ( l2009h presented abundances of Na, a- 
(Si, Ca, and Ti), and iron-peak elements (Fe, and Ni) for a sample 
of seven IC 2391 G/K-type members. They found the follow- 
ing average values: [Fe/H]=-0.01+0.02, [Na/Fe]=-0.02+0.03, 
[Si/Fe]=0.01±0.02, [Ca/Fe]=0.02±0.01, [Ti i/Fe]=0.00±0.02, 
[Ti ii/Fe]=0.I4±0.I3, and [Ni/Fe] =0.00+0.02. Comparing their 
results with our mean abundances (see Table |4) we can conclude 
that the two studies are in excellent agreement; the only exception 
is the Ti abundances derived from the Ti ii lines, for which they 
obtained a 0.12 dex higher value. Note that their [Ti ii/Fe] average 
was based on all sample stars, including those with over-ionisation 
effects. Should we include all sample stars when determining the 
average, then our abundance estimates would agree well. Two of 



iD'Orazi & Randich 


(20091). In Table l5l 


with the ones by 


D'Orazi & Randich 



( l2009l) : the two determinations agree within their respective 
uncertainties for all stellar parameters and final abundances. We 
note that the higher discrepancy for the Ti abundance from the 
Ti II lines is still within the uncertainties, which are significantly 
larger in this case (a different set of line list could partially explain 
su ch a difference ). The b arium abundance for IC 2391 is presented 
in ID'Orazi et al. I ( l2009l) . where they derived a mean abundance 
of [Ba/Fe]=0.68±0.07, which is in excellent agreement with our 
value ([Ba/Fe]=0.62±0.07). 

The chemi cal content of the Argu s association was previously 
determined by IViana Almeida et all ( |2009|) . who gathered Fe, 
Si and Ni abundances for a sample of seven members. They 
also conclude that the association shows a solar abundance 
pattern, finding [Fe/H]=-0.03+0.ofl [Si/Fe]=-0.03±0.03, and 
[N i/Fe] =0.02±0.04. We have four stars in common with the work 
bv lViana Almeida et alj dioO^ . namely CD-28 3434, CD-42 2906, 
TYC 8561-0970-1, and CD-38 583. Difll'erences in Teff are within 
lOOK for all the stars with the exception of TYC 8561-0970-1 
for which we derived a Teff=4900 K, while their final estimate 
is Teff=5348 K. However, ou r final value agrees better with the 
spectral type KO found by [Torres et"^ ( l2006h (see Table [T) 



^ We quote tlieir uncorrected value, and refer the reader to that paper for 
details on the trends with effective temperatures among their sample stai's. 



compared to that by IViana Almeida et al] ( |2009|) . Moreover, 
while gravity values agree well (with an average difference of 
0.125), larger discrepancies are seen for the micro turbulence 
values. IViana Almeida et alj dioO^ derived systematically higher 
microturbulences (on average about 0.15 kms"' higher, which is 
still in fair agreement within the uncertainties) and noted that their 
higher ^ values are probably due to a different set of spectral lines, 
which might be particularly sensitive to the strong magnetic field 
characterising these stars, resulting in larger <f values. Whether 
young pre-main sequence stars exhibit higher microturbulence 
values is still under discussion and there is no general co nsensus 
on this topic jSteenbock & HolwegeJ ( Il98ll) ; |james et alj JToOfii) ). 
Finally, [Fe/H], [Si/Fe] and [Ni/Fe] agree well among the two 
studies, with differences always smaller than 0. 1 dex. 



4 DISCUSSION 



4.1 Chemical tagging 

Stars bom within a single star-forming aggregate share its chemical 
signature with all stars in the cluster, where the elemental abun- 
dance pattern represents the conditions of the proto-cluster gas 
cloud. High resolution studies targeting a large range of chemical 
elements have confi rmed this to be true a mong s e veral Galactic 
open clusters (e.g lOe Silva et ahl l2006l l2007al ; IPancino et al] 
I2OIOI) . Indeed the results presented in Table |4] show that the open 
cluster IC 2391 is also highly chemically homogeneous, with any 
abundance scatter being well within the measurement uncertainty. 

The technique of chemical tagging proposes to use the ele- 
mental abundance patterns of stars to identify dispersed field stars 
with their original formation site dFreeman & Bland-HawthornI 
I2OO2I) . Several recent examples show the concep t of chemical 
taggin g at work; the HR 1 6 14 moving group jPe Silva et alj 
l2007bh . the H ercules stream llBensbv et al .1120071) . the Wolf 630 
moving group ( [Bubar & King' 2010), and th e Hyades Supercluster 
toe Silva et al] [2011; Poinp eia et ah 1201 ih . In all these exam- 
ples the dispersed stellar system was identified via kinematical 
information and its reality as a disrupting cluster was confirmed 
(or disproved in the case of the Hercules stream) via chemical 
abundances. Similarly, the presence of the Argus association was 
identified by the stellar kinematics in the SACY survey and has 
been kinematically associated with the open cluster IC 2391, see 
Figure[T] Probing the chemical information of their stars, we show 
in Table O the average elemental abundances between the open 
cluster stars and the dispersed members of the Argus association, 
are almost identical within the uncertainties. This similarity 
holds for all studied elements, which include those formed via 
the a-process, Fe-peak elements formed via nuclear fusion and 
i-process elements formed via neutron capture. This seems to 
support the case that the stars in the Argus association originally 
formed together with the open cluster IC 2391. Note that unlike 
the previous cases of chemically tagged stellar structures, IC 2391 
and the Argus association are young stellar systems, with an age of 
~30 Myr as presented in Table[3j- 
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Figure 4. [X/Fe] ratios as a function of the effective temperatures, showing the effects of over-ionisation. Note that [Ni/Fe] is not included as the Ni abundances 
do not show any temperature dependent trends (see Table|4)- Circles show IC 2391 stars and triangles show Argus association stars. 
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Table 5. Comparison between our study and the one by D'Orazi & Randich (2009) for two stars in common. 



star 


Teff 


logg 




[Fe/H] 


[Na/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Til/Fe] 


[Till/Fe] 


[Ni/Fe] 




(K) 




(kms-') 
















PMM4362 


5650±50 


4.30±0.10 


1.40±0.15 


-0.05±0.10 


0.09±0.07 


0.05±0.07 


0.11 ±0.09 


0.01±0.11 


0.03±0.11 


0.02±0.11 


PMM4362dro9 


5590±60 


4.45±0.10 


1.15±0.15 


0.00±0.08 


-0.02±0.07 


0.02±0.07 


0.01±0.05 


0.00±0.07 


0.03±0.10 


-0.02±0.07 


PMM4902 


4440±50 


4.50±0.15 


1.30±0.15 


-0.03±0.12 


-0.25±0.15 


0.09±0.08 


-0.25±0.11 


-0.32±0.17 


0.12±0.14 


0.03±0.09 


PMM4902dro9 


4343±60 


4.50±0.10 


1.20±0.15 


0.00±0.09 


-0.34±0.11 


0.06±0.12 


-0.38±0.14 


-0.35±0.17 


0.37±0.14 


0.02±0.12 
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Figure 5. The average abundances for IC 2391 cluster (blue circle) and Argus association stars (red triangle) compared to the solar neighborhood abundances 
of Bensby et al. (2005) (open circles). The average abundance values plotted here do not include the stars which show effects of over-ionisation. Refer to 
section |T4l for details on how the average abundances were derived 



We now compare the cluster and association abundances 
against the solar neighborhood abundances published by 
iBensbv et af] ( |2005|) . As seen in Figure |5] most of the abundances 
of the studied structures sit well within the solar neighborhood 
values. Therefore this lessens the strength of the chemical tagging 
argument discussed above, as other field stars may also share 
a similar chemical abundance pattern as that of the IC 2391. 
Applying the method of chemical t agging we consider the number 
of stars in the iBensbv et all jlOOSh sample which share the same 
abundance pattern as those of IC 2391, with the exception of 
Ba. Keeping a margin of ±0.1 dex we find that 10 stars out of 



the 102 in the IBensbv et al.l ( |2005|) sample could be classed as 
sharing the same abundance patterns (except Ba). If we include 
Ba ab undances in the chem ical tagging process, then only 1 star 
in the IBensbv et alj ( |2005|) sample share a common chemistry as 
IC 2391. Should we adopt a more strict criteria with an abundance 
margin of ±0.05 dex, which is the level of homo g eneity seen within 
bound clusters, then none of the IBensbv et alj ( l2005h stars share 
the open cluster abundance pattern, even without considering Ba. 
Clearly abundance matching of this form is plagued by systematic 
effects between different studies. We await the availability of data 
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from large surveys such a s the GALAH surv ey planned with 
the HERMES instrument JSarden et alj |2010|) . which would be 
homogeneously observed and analysed, to further expand this field 
of research. 



4.2 Ba abundances 

The abundance of Ba, the only analysed .?-process element in 
our study, is clearly over-abundant with respect to the field as 
noted in the discussion above. Similar enhancement in Ba is 
seen in othe r young clusters, including th e Hyades clus t er and 
supercluster dOe Silva et all ilOOd Uoi lh . iD'Orazi et af] ( |2009|) 
detected, in a sample of 20 open clusters, an anti-correlation 
between the [Ba/Fe] ratios and the cluster age. While OCs with 
ages >4 Gyr are characterised by solar Ba abundances, the OCs 
with ages ~ 100-200 Myr exhibit an enhancement up to 0.2-0.3 dex. 
This increasing trend can be reproduced only if the contribution 
from low-mass AGBs (1.0<M/Mo<1.5) to the Galactic chemical 
enrichment is highe r than what w a s pre viously thought (see 
iD'Orazi et all ( |2009|) ; iMaiorca et all feoill) for further details). 
Most intriguingly, the younger clusters of their sample (<70 Myr) 
show an even higher Ba content ([Ba/Fe] up to 0.5-0.6 dex), which 
cannot be explained by this scenario. It is quite unlikely that an 
enrichment in the Ba content took place in the last ~100 Myr 
of Galactic evolution, unless mechanisms o f local enrichment 
are invoked. Similarly, lOesidera et al.l ( l2011a[) . analysing the star 
HD 61005 and its possible link to the Argus association, confirm 
the same abundance pattern, i.e. the [Ba/Fe] is more than a factor 
of four above the solar value. In Figure [6] we report the run of 



[Ba/Fe] with age for clusters present ed in D'O razi et al. (2OO9|0, 



for the three young moving groups bv lD'Or^e t al. (2012), along 
with our measurements for IC 2391 and Argus field stars (this 
study). As one can see, our values fit the Ba-age relationship very 
well, confirming that very young open clusters seem to share an 
extremely high Ba abundance. 

Moreover, the peculiar nature of Ba stands out from the other 
j-process elements: when available, abundances of Y, Zr (first-peak 
j-process element), an d La,Ce (second-peak .^-process element) 
show solar ratios (see iDe Silva et al. I ( l2006l): Ic arrera & Panel nol 
( 1201 ih for the Hvades: ID'Orazi et aljl2oi2l) for three moving 
groups AB Doradus, Carina-Near and Ursa-Major). Since none 
of the current available models can reproduce such a trend 
(accounting for an over-production of Ba without a similar trend 
in La and/or Ce) we are probably dealing with several spurious, 
conspiring effects which mimic a super-solar Ba content. A wide 
discussion of this topic is out of the purpose of the present paper 
(we refer the reader to ID'Orazi et al.l (|2012|) ). here we just recall 
that the presence of hot chromospheres in these young stars might 
play, directly or indirectly, a certain role in the strength of the Ba n 
line(s). 



www. aao. gov. au/HERMES/GAL AH 



4.3 The dissolution of IC 2391 

With all present evidence pointing to a common origin of IC 2391 
and the Argus association stars, we now explor e the likely disso- 
lution mechanism of the original open cluster. iLamers & Gielej 
( l2006l) have shown that most of the energy necessary to dissolve 
a cluster comes from collisions with giant molecular clouds. 
Depending on the impact parameter and the energy involved in the 
collision, a halo of stars is left around a central spheric or more 
oblate core. Qualitatively speaking this is what we see in Figure [T] 
with the IC 239 1 being the core an d the Argus association stars as 
part of the halo dCieles et al.ll200d) . The lack of symmetry around 
the cluster (blue triangles), however, is partially due to the fact tha t 
the IC 2391 is on the limit of the SACY survey dTorres et al.l2008l) 
of 150 pc. The second half of the halo being beyond 150 pc. 

For a small-N system like IC 2391 a halo of unbound stars 
can also be the result of ejections of stars from the cluster core 
due to two-body encounters. Even when the tidal field is weak, 
clusters lose stars because of close encounters of stars in the 
cluster. This process is important for clusters whose age exceeds 
the half-mass relaxatio n time Tfh. Based on an estimated half-mass 
radius of about 1.5 pc i Dias et al]|2002l) . a mass of about 500 M© 
dBoudreault & Bailer- Jonesll2009l) we use t he standard expression 
for the half-mass relaxation time-scale Tph of lSDitzer&Har3dl97ll) 
and find T,h - 70 Myr. Here we assumed the mean mass of the 
stars to be 0.5 Mq (i.e. A' = 1000) and we use 0.02A' as the 
argument of the Coulomb logarithm. Although this value for T|h is 
uncertain to within a factor of two or three, the similarity between 
T,!, and the age of IC 2391(about 30+10 Myr) suggests that it is 
possible that IC 2391 is dynamically evolved and has ejected a 
fraction of its stars. These unbound stars we now observe as the 
Argus association around the c luster. It implies that the cluster was 
smaller and denser in the past dHenonll 19651 : iGieles et al.ll2010l) . 

In fact, we can speculatively estimate the initial size and 
velocity dispersion of the cluster. The radius of the Argus asso- 
ciations is about 150 pc. Combined with an age of 30 Myr for 
IC 2391 this means that the most distant stars must have escaped 
with a velocity of about 5 kms"'. The typical velocity with which 
stars escape from an isolated cluster is few Vi^s, where v,n,s is 
the root-mean square velocity of cluster. So the most distant 
stars escaped when the cluster had i^js ^1-2 kms"'. Combined 
with the mass this corresponds to a half-mass radius of about 
Th-0. 2-0.4 pc. In Fig.jTjwe illustrate this idea with the result of an 
A'-body simulation of an isolated cluster with an initial rh=0.25 pc. 
After 30 Myr there are unbound stars out to 100 pc away from 
the cluster, which have velocities of a several kms"'. The current 
ri, 1 . 1 pc and 897 stars are still bound to the cluster, comparable 
to IC 2391. 

This implies that the initial density of IC 2391 was several 
orders of magnitude higher than the present day density and the 
initial relaxation time a factor of a few shorter. This could be 
the explanation for the observed mass segregation in the cluster 
dBoudreault & Bailer-Jon"esll2009l) . These are of course very rough 
estimates, because we have ignored the presence of the Galactic 
tidal field and external factors such as giant molecular clouds. 



^ Due to a systematic offset between dwarfs and giants in the lD'Orazi et al.l 
)2009h sample, here we consider, for homogeneity, only clusters whose 
abundances come from un-evolved members. 
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Figure 6. [Ba/Fe] as a function of cluster age for IC 2391 (filled circle) and the Argus association (filled triangle) as well as open clusters (labelled as open 
squares from D'Orazi et al. 2010) and moving groups (labelled as open circles from D'Orazi et al. 2012). 



5 SUMMARY AND CONCLUSIONS 

In this paper, we presented the members of the Argus association 
as identified via the SACY survey and find the age of the system 
to be about 30±10 Myr. Based on high resolution UVES spectra 
we determined the abundances for Fe, Na, Mg, Al, Si, Ca, Ti, 
Cr, Ni and Ba in both the open cluster IC 2391 and the Argus 
association. All stars in the open cluster and Argus association 
were found to share similar abundances with the scatter well within 
the uncertainties, where [Fe/H] = -0.04 ± 0.03 for cluster stars and 
[Fe/H] = -0.06 ± 0.05 for Argus stars. All other elements were in 
their solar proportions with the exception of Ba which we observed 
to be enhanced to around 0.6 dex. We also discussed the effects of 
over-ionisation/excitation observed in our sample for stars cooler 
than about 5200 K, as has been previously noted in the literature. 

In summary we present strong kinematic, evolutionary and 
chemical evidence to support a common origin for the Argus 
association and the open cluster IC 2391. In particular the chemical 
tagging of the stars presented in this paper provides a necessary 
condition to prove the link between the Argus association and 
IC 2391. However, as discussed in the text, many elemental 
abundances of the studied stars are indistinguishable from the solar 
neighborhood, with the exception of Ba which is an important 
piece of evidence in favor of a common origin of the two stellar 
systems. For reasons which are not yet completely understood, 
barium abundances have been recently found to correlate with 



age. The abundances found for both stellar systems are similar. 
Moreover, they fit the Ba-age correlation very well. 

Given all evidence in hand, we conclude that we are witness- 
ing the on-going dissolution of IC 2391, where the stellar mem- 
bers of the Argus association stars were originally born in the 
same proto-cluster cloud as IC 239 1 . Simple modeling of this sys- 
tem finds the dissolution of this timescale to be consistent with 
two-body interactions. A more detailed modeling of the disruption 
mechanisms, predicting spatial distribution and possibly the veloc- 
ity dispersion of both groups would strengthen the proposed sce- 
nario. From an observational point of view, the second half of the 
stellar halo will be easily found by GAIA settling the question. 
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Figure 7. Results of an N-body simulation with NBODY6 (Aarseth 2003) of an isolated cluster with N=1024, a Kroupa (201 1) stellar mass function between 
0.1 Mq and 10 Mq. The initial density profile followed from a Plummer (1911) model with a half-mass radius of r), = 0.25 pc. Small blue symbols represent 
all bound stars at an age of 30 Myr and the larger red circles are unbound members. The cluster and halo are comparable to IC 2391 and the Argus association 
at 30 Myr. 
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